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Abstract. 

A search for the double beta decay of zinc and tungsten isotopes has been performed 
with the help of radiopure ZnW04 crystal scintillators (0.1 — 0.7 kg) at the Gran 
Sasso National Laboratories of the INFN. The total exposure of the low background 
measurements is 0.529 kg x yr. New improved half-life limits on the double beta 
decay modes of ^''Zn, ^"Zn, ^^°W, and ^^^W have been established at the level of 
10^^ — 10^^ yr. In particular, limits on double electron capture and electron capture with 
positron emission in ^^Zn have been set: T'lj^'^ > 1.1 x lO^^ yr, T°^^^ > 3.2 x 10^° yr, 

T^jf^ > 9.4 X 10^0 yr, and T°jf^ > 8.5 x 10^" yr, ah at 90% C.L. Resonant 
neutrinoless double electron capture in ^^"W has been restricted on the level of 
T°j2^ > 1.3 X 10^^ yr. A new half-life limit on a transition of ^^-^W to the metastable 
excited level 1/2" 375 keV of "^Hf has been estabhshed: T1/2 > 6.7 x 10^" yr. 
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1. Introduction 

Double beta (2/3) processes are nuclear transformations when the charge of nuclei 
changes by two units: [A, Z) — {A, Z ±2). There are two main modes of 2^ 
decay: two neutrino mode (2z/) when two neutrinos are emitted together with two 
beta particles, and neutrinoless mode (Oz/). Oz/2/3 decay violates the lepton number by 
two units and therefore is forbidden in the Standard Model (SM) [Ij. However, the 
Oi^2/3 decay is predicted in some SM extensions where neutrino is expected to be a 
true neutral particle equivalent to its antiparticle (Majorana particle) [2]. Experiments 
on neutrino oscillations already gave evidence for neutrino to be massive |3], however 
these experiments are sensitive only to the differences of squared masses of neutrinos. 
The observation of Oz/2/3 decay could resolve important problems of particle physics: 
what is the absolute scale of neutrino mass? Which neutrino mass hierarchy (normal, 
inverted, or quasi-degenerate) is realized in nature? Is the neutrino Majorana (z/ = v) 
or Dirac (z/ ^ v) particle? Is the lepton number absolutely conserved? Additionally, 
investigations of neutrinoless double (3 decay could test admixture of right-handed 
currents in electroweak interaction and existence of majoron^. 

While 21/2(3 decay is allowed in the SM, it is a second order process in perturbation 
theory characterized by extremely low probability. Investigations of the 2v2(3 decay 
examine theoretical calculations of the nuclear matrix elements, contributing to the 
development of theoretical description of Oz/2/3 decay. 

Double beta decay experiments are concentrated mainly on 2(3 processes with 
emission of two electrons (2/3~). Two neutrino mode of 2(3~ decay was detected for 11 
nuclides among 35 candidates; corresponding half-lives are in the range of 10^® — 10^^ yr 
P Uni HB 112]. In addition, the 2^2(3' transitions of ^°°Mo and ^^°Nd to the first 0+ 
excited states of daughter nuclei were observed too [U [101 E]- To-date the 2z/2/3~ 
decay is the rarest radioactive decay ever discovered. Developments in the experimental 
techniques during last two decades lead to impressive improvement of sensitivity to the 
neutrinoless mode of 2(3~ decay up to the level of T1/2 ~ 10^^ — 10^^ yr P [TT] . Moreover, 
some possible positive indication for ^^Ge with T^l2^ = 2.2 x 10^^ yr has been mentioned 
in [13], and new experiments are in preparation both on ^^Ge [HI [15] and other isotopes. 

A more modest sensitivity was reached in the experiments searching for 2/3 processes 
with decreasing charge of nuclei: capture of two electrons from atomic shells (26), 
electron capture with positron emission (e/3+), and double positron decay (2/3+). There 
are 34 possible candidates for 2e capture; among them, only 22 and 6 nuclei can also 
decay through e(3~^ and 2/3"'' channels, respectively [9]. In contrast to the 2/3~ decay, 
even the allowed 2z/ mode of 2e, e(3^, and 2/3+ processes are still not detected in 
direct experiment^ and the obtained half-life limits are much more modest. The most 

§ Massless or light bosons that arise due to a global breakdown of {B — L) symmetry, where B and 
L are the baryon and the lepton number, respectively. Literature considers Oi^2/3 decay channels with 
one (0i^2^Ml) I1[S], two {Qv2(3M2) 017], and "bulk" {Qv2(3hM) [5] majoron emissions. 
II For completeness, we remind that a possible evidence of 2v2£ capture in ^^°Ba with T"^^^^ ~ 
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sensitive experiments have given limits on the 2e, e(3^ , and 2/3+ processes at the level 
of 10^^ — 10^^ yr [9*, TT]. Reasons for such a situation are: 1) lower energy releases 
in comparison with those in 2(3~ decay, that results in lower probabilities of the 
processed, as well as provides difficulties to suppress background; 2) usually lower 
natural abundances (5) of 2/3"*" isotopes (which are typically lower than 5% with only 
a few exception^). Nevertheless, studies of 2e and el3^ decays are important, because 
observation of neutrinoless mode of such process could help to distinguish between the 
mechanisms of neutrinoless 2/3 decay (is it due to non-zero neutrino mass or to the 
right-handed admixtures in weak interactions) |18j . 

Zinc tungstate (ZnW04) scintillators contain four potentially 2/3 active isotopes: 
e^Zn, ™Zn, ^«°W, and ^^^W (see Table 1). It is worth to mention, ^^Zn and ^^^W 
have comparatively large natural abundance that allows to apply ZnW04 detectors 
without high cost enriched isotopes. Moreover, the 2z/2/3~ decay of ^^^W is expected 
to be strongly suppressed [19], that could provide favorable conditions to search for 
neutrinoless 2/3~ decays, including processes with emission of majoron(s) which have 
broad energy spectra, somewhat similar to that of the two neutrino mode. The ^^°W 
isotope is also an interesting 2/3 nuclide because, in case of the capture of two electrons 
from the K shell {E^ = 65.4 keV), the decay energy is rather small (13 ± 4) keV. Such 
a coincidence could give a resonant enhancement of the double electron capture to 
the corresponding level of the daughter nucleus [201I2T1 [22| [23t [2il [2it [25] . 

Table 1. Potentially 2/3 active isotopes of zinc and tungsten present in ZnW04 crystal 
scintillators. 



Transition 


Energy release 
{Qpp), keV m 


Isotopic 
abundance 

(%) m 


Decay 
channels 


Number of mother 
nuclei in 100 g of 
ZnW04 crystal 




1095.7(0.7) 


49.17(75) 


2£, e/3+ 


9.45 X 10^2 


™Zn -> ^OGe 


998.5(2.2) 


0.61(10) 


2/3- 


1.17 X 10^1 


180^ _^ ISOjjf 


144(4) 


0.12(1) 


2e 


2.31 X 10^0 


186^ ^ 1860s 


489.9(1.4) 


28.43(19) 


2/3- 


5.47 X 10^2 



The best to-date half-life limits on different modes and channels of 2(3 processes in 
zinc and tungsten isotopes (except of Oz/2/3- decays of ^*^W) were obtained in previous 
stages of this experiment [281 129]. The best half-life limits on Oz/2/3- decays of ^^^W 
to the ground and excited states of ^^^Os were set in the Solotvina experiment with 
cadmium tungstate scintillator enriched in ^^^Cd [30] . 

(0.5 — 2.7) X 10^^ yr has been reported in geochemical studies [TBI [TT] . 

% The value of half- life is inversely related to the phase-space factor (G); the latter depends on energy 
release as G '-^ Q^^j for 2i^2/3 decay and ~ for Qv2(3 decay 

+ Only 6 nuclides from a complete list of 34 isotopes-candidates on 2e, £(i^ , and 2/3+ processes have 
natural abundances of more than 5% [5]. 
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Here we present the final results of the experiment to search for double beta 
processes in zinc and tungsten with the help of ZnW04 crystal scintillators. As a 
by-product of the experiment, we also have set a new limit on a decay of ^^•^W to the 
375 keV metastable excited level of ^^^Hf. 

2. Experiment and data analysis 

The low background experiments to search for double beta processes in zinc and 
tungsten isotopes have been performed by using zinc tungstate crystal scintillators. The 
scintillation detectors with ZnW04 crystals, the experimental set-up, the measurements 
and the data analysis are described in detail in [28l |29l [31] . Here we outline the main 
features of the experiment. 

2.1. ZnW04 crystal scintillators 

Four ZnW04 crystal scintillators were used in our studies. Two crystals (117 g, 
20 X 19 X 40 mm, and 699 g, 044 x 55 mm) were produced by the Czochralski method 
[321 133] ill the Institute for Scintillation Materials (Kharkiv, Ukraine). After 2130 h of 
low-background measurements the crystal of 699 g was re-crystallized with the aim to 
study the effect of the re-crystallization on the radioactive contamination of the material. 
The third ZnW04 crystal (141 g, 027 x 33 mm, the sample had slightly irregular shape) 
was obtained by the re-crystallization process and used in further measurements. The 
fourth ZnW04 crystal scintillator (239 g, 041 x 27 mm) was produced in the Nikolaev 
Institute of Inorganic Chemistry (Novosibirsk, Russia) by the low-thermal gradient 
Czochralski technique [Ml [35] . The radioactive contaminations of the used crystals 
are reported in [31] . 

2.2. Low-background measurements 

The ZnW04 crystal scintillators were fixed inside a cavity of 049 x 59 mm in the central 
part of a cylindrical polystyrene light-guide of 066 x 312 mm. The cavity was filled up 
with high purity silicone oil. The light-guide was optically connected on opposite sides by 
optical couplant to two low radioactivity EMI9265-B53/FL 3" photomultipliers (PMT). 
The light-guide was wrapped by PTFE reflection tape. The detector was modified at 
the final stages of the experiment: two polished quartz light-guides (066 x 100 mm) 
were installed between the polystyrene light-guide and the PMTs to suppress 7 ray 
background from the PMTs. 

The detector has been installed in the low background DAMA/R&D set-up at the 
underground Gran Sasso National Laboratories of the INFN (Italy) at the depth of ~ 
3600 m w.e. It was surrounded by Cu bricks and sealed in a low radioactive air-tight 
Cu box continuously flushed with high purity nitrogen gas (stored deep underground 
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for a long time) to avoid presence of residual environmental radon. The outer passive 
shield consisted of 10 cm of high purity Cu, 15 cm of low radioactive Boliden lead, 1.5 
mm of cadmium and 4/10 cm polyethylene/paraffin to reduce the external background. 
The whole shield has been closed inside a Plexiglas box, also continuously flushed by 
high purity nitrogen gas. An event-by-event data acquisition system accumulates the 
amplitude, the arrival time, and the pulse shape of the events. 

The energy scale and the energy resolution of the ZnW04 detectors have been 
measured with 7 sources ^^Na, ^°Co, -'^'^^Ba, -'^'^'^Cs, ^^^Th, and ^'^-'^Am. The energy 
resolution of the detectors (full width at half of maximum) was in the range of 
(8.8 - 14.6)% for 662 keV 7 line of ^^^Cs. 

2. 3. Interpretation of the background 

As an example, the energy spectrum accumulated over 4305 h with the 239 g ZnW04 
crystal scintillator in the low background set-up is shown in Fig. 1. The energy spectrum 
accumulated over 2798 h in the same conditions with low (~ 10 keV) energy threshold 
is presented in Inset. A few visible peaks in the spectrum can be ascribed to 7 quanta 
of naturally occurring radionuclides "^^K, ^^^Bi (^^^U chain) and ^os^jj (^^^Th) from the 
materials of the set-up. The presence of the peak with energy ^ 50 keV can be explained 
by internal contamination of the crystal by ^^°Pb. A comparatively wide peculiarity at 
the energy 0.8 MeV is mainly due to the a active nuclides of U and Th chains present 
in the crystal as trace contamination. 

The radiopurity of the ZnW04 scintillators was already estimated [291 EI] by using 
the data of the low-background measurements. The time-amplitude analysis (see details 
in [MIET]), the pulse-shape discrimination between /9(7) and a particles |38j, the pulse- 
shape analysis of the double pulses (overlapped Bi-Po events) [HOI |39l |40] , ^-iid the Monte 
Carlo simulation of the measured energy spectra were used to determine radioactive 
contamination of the ZnW04 crystals. The radioactive contamination of the ZnW04 
crystals is on the level of 0.002 — 0.8 mBq/kg (depending on the source); the total a 
activity is in the range 0.2 — 2 mBq/kg. Moreover, particular contaminations associated 
with the composition of ZnW04 detector were observed [31]: the EC active cosmogenic 
(or/and created by neutrons) nuclide ^^Zn (T1/2 = 244.26 d [H]) with activity 0.5 — 0.8 
mBq/kg (depending on the ZnW04 sample) and the a active tungsten isotope ^*°W 
(with half-life: T1/2 ^ 10^^ yr [3ll [391 ill Il3] , and energy of the decay: Q„ = 2508(4) keV 
[26] ) with activity 0.04 mBq/kg (see Fig. 1). 

3. Results and discussion 

3.1. Response of the ZnWO/^ detectors to 2(3 processes in zinc and tungsten 

The response functions of the ZnW04 detectors for the 2/3 processes in Zn and W 
isotopes were simulated with the help of the GEANT4 package [H] with the Low Energy 
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Figure 1. (Color online) The energy spectrum accumulated with the ZnW04 crystal 
scintillator 041 x 27 mm in the low background DAMA/R&D set-up over 4305 h. The 
energy spectrum of a events selected by the pulse-shape discrimination is drawn by 
points. Fit of the a peak of ^^"W by Gaussian function (solid line) is shown. (Inset) 
The energy spectrum of 7 and j3 events selected by the pulse-shape discrimination 
technique from the data measured over 2798 h with the same crystal scintillator in the 
set-up with lower energy threshold and with additional quartz light-guides. Energies 
of 7 lines are in keV. 

Electromagnetic extensions. The initial kinematics of the particles emitted in the decays 
was generated with the DECAYO event generator [15]. As examples, the expected energy 
distributions for the ZnW04 detector 044 x 55 mm are shown in Fig. 2 and Fig. 3. The 
background models included the internal contamination of the ZnW04 scintillators (^*^K, 
^°Co, ^^Zn, ^''Rb, ^'^Sr-^'^Y, ^^'^Cs, active nuclides from U/Th families), and the external 
7 rays from radioactive contamination of the PMTs and the copper box (^''K, ^^^Th, 
^^^U); they were also simulated with the help of the GEANT4 and DECAYO packages. 

3.2. Double (3 processes in ^^'™Zn and i80,i86 ^ 

Comparing the simulated response functions with the measured energy spectra of 
the ZnW04 detectors, we have not found clear peculiarities, which can be evidently 
attributed to double beta decay of zinc or tungsten isotopes. Therefore only lower 
half-life limits can be set according to the formula: 



lim Ti/2 = N ■ rj ■ t ■ ln2/ lim S, 



(1) 
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Figure 2. (Color online) Simulated response functions of the detector based on the 
ZnW04 scintillator 044 x 55 mm for the different 2/3 processes in Zn isotopes: (a) 2e 
capture in ^^Zn; (b) e/S"*" decay of ^'^Zn; (c) 2/3~ decay of '''^Zn; (d) Qu2(3~ processes 
with majorons emissions in "^^Zn. One million decays was simulated for each process. 



where N is the number of potentially 2/3 unstable nuclei in a crystal scintillator, rj is 
the detection efficiency, t is the measuring time, and limS* is the number of events of 
the effect searched for which can be excluded at a given confidence level (C.L.; all the 
limits in the present study are given at 90% C.L.). 

For the 2z/ double electron capture in ^^Zn from the K shell, the total energy released 
in the detector is equal to 2Ek = 16.7 keV (where Ek is the binding energy of electrons 
on the K shell of nickel atoms). The detection of such a small energy deposit requires 
rather low energy threshold. In our measurements with the ZnW04 crystal scintillator 
041 X 27 mm the energy threshold of 10 keV was enough low (see Fig. 1, Inset) to 
observe at least the higher energy part of the 2v2K peak. Moreover, the background 
level (which is mainly due to PMT noise in the low energy region) was decreased in 
comparison to our ffist measurement [28] thanks to the improved scintillation properties 
of the ZnW04 crystal (slightly higher transmittance, light output and energy resolution) 
and the enhanced light collection from the scintillator. The light collection was increased 
by special treatment of the crystal surface, which was diffused with the help of grinding 
paper (in our first experiment, the ZnW04 crystal scintillator was polished |28]). Finally, 
a significant difference of ZnW04 pulse-shape (effective decay time is ~ 24 /is ^6]) in 
comparison to much faster PMT noise (few nanoseconds) offers the possibility to exploit 
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Figure 3. (Color online) Simulated response functions of the ZnW04 detector 
044 X 55 mm for the different 2/3 processes in W isotopes: (a) 2e capture in ^^^W; (b) 
and (c) 2/3~ decay of -'^^^W to the ground and excited states of ^^^Os, respectively; (d) 
Oi/2/3~ decays of -'^^^W with majorons emissions. One million decays was simulated for 
each process. 



the rejection of residual PMT noise by using the pulse-shape discrimination. However, 
this procedure eliminates some part of scintillation signals near energy threshold. The 
energy dependence of the detection efficiency was determined with the help of ^^^Ba, 
^^''Cs, ^^^Th, and ^^^Am radioactive sources. The measured efficiency ranges from about 
55% at 15 keV up to about 95% at 30 keV (one can compare these values with the 
detection efficiencies 30% at 15 keV and 65% at 30 keV obtained in p8]). 

To set a limit on the 2v2K decay of ^"^Zn, taking into account the proximity of the 
energy threshold and the contribution from remaining PMT noise, we use a conservative 
requirement: the theoretical energy distribution should not exceed the experimental one 
in any energy interval, including error bars in the experimental values (see Fig. 4). In 
this way the limit on the peak area is lim S = 4665 counts. Taking this value (already 
corrected for the efficiency) for the peak area, we conservatively give the following half- 
life limit on the 2u2K process: 

T2-2^(64Zn) > 1.1 X W yr. 

To estimate limits on other double /3 processes, we have used the following approach: 
the energy spectrum was fitted in the energy range of an expected 2/3 signal by a model 
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Figure 4. The energy spectrum of the ZnW04 crystal scintillator 041 x 27 mm 
measured over 2798 h, corrected for the energy dependence of detection efficiency, 
together with the 21^2^ peak of ^^Zn with T'lj^'^ = 1.1 x yr excluded at 90% C.L. 



built by the simulated distributions of internal and external background and of the effect 
searched for. The background model was composed of '^''K, ^^Zn, ^''Sr-^^Y, ^^"^Cs, U/Th 
inside a crystal (for fit of a low energy part of the data we have also used a model of 
internal ^^Rb), and "^^K, ^'^^Th, ^ss^j ^j^g PMTs and the copper box. The activities of 
the U/Th daughters in the crystals have been restricted taking into account the data on 
the radioactive contamination of the ZnW04 crystal scintillators [HT] . The initial values 
of the ^°K, ^^^Th and ^^'^U activities inside the PMTs have been taken from Ref. 
while activities inside the copper box have been assumed to be equal to the estimations 
obtained in Ref. [H]. We have used different combinations of the accumulated data to 
reach the maximal sensitivity to the double beta processes searched for. Additionally 
we have also applied the so called la approach when a statistical uncertainty of the 
number of events accumulated in the energy region of the expected 2(3 signal (square 
root of the number of events) was taken as lim S. This simple method allows to obtain 
a correct evaluation of the experimental sensitivity to the 2(3 process searched for. It 
should be stressed that the detection efficiencies in all the distributions analyzed are 
at least 99.9% for all the processes. Taking into account the efficiency of '^{(3) events 
selection by the pulse-shape discrimination (98%), the total detection efficiencies are at 
least 97.9% for all the 2(3 processes searched for. 

Let us give an example of the analysis by using the two approaches to search for 
electron capture with positron emission in ^^Zn. 14922 events were observed in the 
energy interval 530 — 1190 keV of the spectrum accumulated with an exposure 0.3487 
kg X yr (see Fig. 5), which gives limS* = 122 counts. With the detection efficiency 
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in the energy interval to the 2veP^ decay of ^^Zn (82%), one obtains the half-hfe hmit 
Ty2^ > 1.5 X 10^^ yr at 68% C.L. In order to apply the second approach, the starting 
and final energies of the fit were varied as 380—550 keV and 1260—1430 keV, respectively, 
with the step of 10 keV. The result of the fit in the energy region 520 — 1350 keV was 
chosen as final giving the minimal value of l^-d-f- = 119/98 = 1.23. It gives the total 
area of the 2vef5^ effect (—208 ±254) counts which corresponds (in accordance with the 
Feldman- Cousins procedure [19]) to \im.S = 238 counts in the full energy distributions 
for 2z/£:/3+ decay. Thus, one can calculate the following half-life limit, rather similar to 
the value obtained by using the la approach: 




500 1000 1500 201 

Energy (keV) 



Figure 5. (Color online) The measured energy spectrum of the ZnW04 scintillation 
crystals (the total exposure is 0.349 kg x yr) together with the GEANT4-simulated 
response functions for e/3+ process in ^"^Zn excluded at 90% C.L. The most important 
components of the background are shown. The energies of 7 lines are in keV. 



Tyf(64Zn) > 9.4 X 1020 yr. 

In case of the neutrinoless electron capture with positron emission, the spectrum 
with the total exposure 0.3487 kg x yr was fitted in the energy interval (410 — 1370) 
keV {x^/n.d.f. = 113/94 = 1.2). The fit gives the area of the effect searched for as 
(52 ± 129) counts, which corresponds (again in accordance with the Feldman-Cousins 
procedure) to lim S = 264 events. It allows to set the following limit on the half-life of 
^^Zn relatively to the 0i>e(3^ decay: 

T°;f^(64Zn) > 8.5 X 1020 yr. 

The energy distributions expected for the 2i/e(3^ and Oz/£/3+ decay of ^^Zn, excluded at 
90% C.L., are shown in Fig. 5. 
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In case of 0z/2e decay of ^"^Zn, different particles are emitted: X rays and Auger 
electrons from deexcitations in atomic shells, 7 quanta and/or conversion electrons from 
deexcitation of daughter nucleus. We suppose here that only one 7 quantum is emitted 
in the nuclear deexcitation process; it is the most pessimistic scenario from the point 
of view of registration of such an event in a peak of full absorption at the Qpp energy. 
Unfortunately, 2K, KL, 2L (and other) modes are not energetically resolved in the high 
energy region due to finite energy resolution of the ZnW04 detectors. So, the fit of the 
measured spectrum (exposure 0.3647 kg x yr) in the energy interval 440 — 1350 keV 
{x^/n.d.f. = 98/89 = 1.1) gives the area of the 0^26 effect searched for as (—780 ±853) 
counts. Taking into account the Feldman-Cousins procedure, we calculated lim 5* = 742 
events and the following limit on 0^26 transition of ^^Zn to ground state of ^^Ni: 

T0/|^(64Zn) > 3.2 X 10^° yr. 

Limits on double electron capture in ^^°W were set by analyzing all the data 
accumulated in the experiment over 0.529 kg x yr. The low energy part of the spectrum 
is shown in Fig. 6. The least squares fit of this spectrum in the 100 — 260 keV energy 
interval gives (141 ±430) counts for the 2v2K peak searched for {x^/n.d.f. = 5.39/5 = 
1.08), providing no evidence for the effect. These numbers lead to an upper limit of 846 
counts. Taking into account the detection efficiency for this process close to 98%, one 
can calculate the half-life limit: 

r2;|^(i8ow) > 1.0 X 10^8 yr. 

The same approach gives the limit for the neutrinoless 2e process in ^^''W: 

^0.2.(180^^) > 1.3 X 10^8 

The expected energy distributions for 0z/2e and 2v2K decay of ^^°W corresponding 
to the best previous restrictions obtained in the Solotvina experiment [30] with the 
help of low background cadmium tungstate crystal scintillators are presented in Fig. 6. 
Advancement of the sensitivity in the present study was reached thanks to the lower 
background of ZnW04 detectors in comparison to CdW04 where the counting rate in 
the energy interval up to 0.4 MeV was caused mainly by the /3 decay of ^^^Cd. The 
Qu2e decay of ^'^°W is of particular interest due to the possibility of resonant process 

By using the approaches described above, the half-life limits on other 2(3 decay 
processes in ^^Zn, ^''Zn, and ^^^W have been obtained. All the results are summarized 
in Table 2, where the data of the most sensitive previous experimental investigations 
and theoretical estimations are given for comparison. 

The obtained bounds are well below the existing theoretical predictions; 
nevertheless most of the limits are higher than those established in previous experiments. 
It should be stressed that in contrast to the results obtained in researches of double /3~ 
decay (sensitivity of the best experiments is on the level of 10^^ — 10^^ yr [9l [TOl [TT]). 
only five nuchdes (^°Ca [50], ^^Zn [present work], ^^Kr [51], ^^^Sn [52], and ^^^Te [53]) 
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Figure 6. (Color online) Energy spectrum of the background of the ZnW04 detectors 
(exposure 0.529 kg x yr). The simulated response functions for double electron capture 
in W are shown; the half-lives Tfj^'^ = IxW^ yr and T^J^'' = 9x IQi^ yr correspond 
to the best previous limits obtained in [30] with the help of cadmium tungstate crystal 
scintillators. 



among 34 potentially 2e, and 2/3+ active isotopes were investigated at the level of 
sensitivity limTi/2 ~ 10^^ yr. 

3.3. Search for a decay of tungsten isotopes 

In addition to the previous observation of the a decay ^'^^W -> ^^^Hf (g.s. to g.s. 
transition) with CdW04 and CaW04 detectors [391 1321 SS], this rare process was 
observed also in our data with ZnW04 scintillators with T1/2 = 1.310:5 x 10^^ yr [31] 
(one can also see the a peak of ^'^''W in the a spectrum presented in Fig. 1). 

Here we report a new limit on the a decay of ^^'^W {Q^ = 1680(2) keV [2S]; 5 = 
14.31(4)% [27]) to the 1/2" metastable level of ^^^Hf (375 keV, T1/2 = 18.67 s [H]). 
The search for this process has been performed by using the data of all the runs with 
the ZnW04 detectors with the total exposure 0.5295 kg x yr. The signature of such 
a transition is delayed 7 quanta after the emission of the a particle. The expected 
distribution of the time intervals between the a and the 7 events should correspond to 
Ti/2 = 18.67 s. The time-amplitude technique [361 137] and the pulse-shape discrimination 
method [381 HH] have been applied to search for the a decay. Taking into account the 
a/l3 ratiqj {a/(3 ^ 0.17) for the ZnW04 scintillator [46j, we expect to observe the a 
peak of the ^'^^W decay to the ^^^Hf metastable level at the energy 220 keV in 7 scale, 

* It is defined as ratio of a peak position in the energy scale measured with 7 sources to the real energy 
of a particles. 



Final results of experiment to search for 2(3 processes in zinc and tungsten ... 13 



Table 2. Half-life limits on 2/3 processes in Zn and W isotopes and comparison with 
the theoretical predictions. Quoting best previous experimental results, we exclude 
limits obtained on previous stages of our experiment [551 HI] • 



Transition 


Decay 
channel 


Level 

of daughter 
nucleus 


Experimental limits on T1/2, yr at 90% C.L. 


Theoretical estimations of 

the half-lives T1/2, yr 

((■m„) = 1 eV for 0^2/3 decay) 


Present work 


The best previous results 


"■'Zn ^ '■'■'Ni 


2u2K 
Qu2e 


g-s. 
g-s. 


> 1.1 X IQi" 

> 3.2 X 10^° 


> 6.0 X 10^'' H 

> 7.4 X lO'** 


(1.9-7.1) xio™ mi 
(1.2 ±0.2) X inzi 


2uel3+ 


g-s. 
g-s. 


> 9.4 X 10™ 

> 8.-5 X 10™ 


= (1.1 ±0.9) xioi" Eg 

> 1.3 X 10™ ISHI 

> 1.3 X 10™ mi 


(0.9-2.2) xlO-''' 
(4.7 ±0.9) x 10'" |S7| 


™Zn ^ ™Ge 


2u2p- 
Qu2p- 

au2p-Mi 

Qu2l3-M2 
0u2p-bM 


g-s. 

g-s. 
g-s. 
g-s. 
g-s. 


> 3.8 X 10'** 

> 3.2 X 10'" 

> 6.0 X 10'*^ 

> 4.7 X 10'** 

> 5.4 X 10'** 


> 1.3 X lo'*' im 

> 7.0 X 10'^ uni 


4.5 X 10" - 3.6 X 10^' mi 
2.5 X 10" - 6.4 X 10^=* inn 

7.0 X 10™ mi 
> 3.1 X 10^^ inii 
9.8 X 10™ mi 




2v2K 
{)u2s 


g-s. 
g-s. 


> 1.0 X 10'** 

> 1.3 X 10'** 


> 7.0 X lo'^' im 

> 9.0 X lo'" mi 


2.5 X 10^' - 2.5 X 10™ [23 
3.3 X 10^" - 5.0 X lO'™ P4 
3.0 X 10™ - 4.0 X 10^^ psj 


180 w ^ is^Os 


2v2B- 

2u2p- 

Qu2l3- 

Qu2l3- 

0u2p-Ml 

Qu2p-M2 

Qv2B-hM 


g-s. 

2+ (137 keV) 
g-s. 

2+ (137 keV) 

g-s. 

g-s. 

g-s. 


> 2.3 X lO'" 

> 1.8 X 10™ 

> 1.0 X 10^1 

> 9.0 X 10™ 

> 5.8 X lO'" 

> 1.1 X lO'" 

> 1.1 X lO'" 


> 3.7 X lo'** mi 

> 1.0 X lo'" mi 

> 1.1 X 10^1 mi 

> 1.1 X 10^1 mi 

> 1.2 X 10™ mi 


7.1 X 10™ - 1.2 X 10™ Pi 
> 6.1 X 10^* 113 

6.4 X 10^" mi 



with energy resolution FWHMq, = 62 keV. All the a events selected within 150 — 270 keV 
have been used as triggers, while a time interval 0.1 — 60 s (88.9% of ^''^Hf* decays) and 
a 325 — 425 keV energy window have been set for the second 7 events (energy resolution 
for gammas at the energy 375 keV: FWHM^ = 64 keV). Ninety five pairs were selected 
from the data. The fit of the distribution of the selected "a events" by a simple model 
built by a first degree polynomial function (to describe the background) plus a Gaussian 
(the a peak searched for) gives the area of the effect searched for as (10.5 ±17.6) counts, 
which corresponds to limS* = 39.4 events. The excited 375 keV level of ^^^Hf deexcites 
with emission of two 7 quanta of 161 keV and 214 keV [H]. The efficiency to detect a 
peak at the total energy release of 375 keV in ZnW04 detectors was calculated with the 
GEANT4 |13]; it was equal 0.71 to 0.86 in dependence on the volume of the ZnW04 
crystal. The half-life limit was calculated according to the formula analogous to (1): 

lim Ti/2 = In 2 ■ ripsn ■ T.^! ■ N ■ t/ YrniS , 

where rjpsD is the efficiency of the pulse-shape discrimination (37.2%), is the number 
of ^^^W nuclei, rj is the registration efficiency of the total energy release of 375 keV, 
and t is the time of measurements with specific ZnW04 detector. In result, we set the 
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following limit on the half-life of the a decay of ^^^W to the metastable 375 keV excited 
level of i™Hf: 

Ti/2 (^^^W ^ i^^Hf*, 375 keV) > 6.7 x lO^o yr. 




Energy (keV) 

Figure 7. (Color online) Energy spectrum of the events selected by the time-amplitude 
and the pulse-shape analyses from the data accumulated by ZnW04 detectors with an 
exposure 0.5295 kg x yr. These events satisfy the search criteria for a transition of 
183^ to metastable level of ^^^Hf. Polynomial function used as a background model 
and the Gaussian peak corresponding to a decay of ^^■^W with the half- life T1/2 = 
6.7 X 10^° yr excluded at 90% C.L. are also shown. 

The energy spectrum of the selected events is shown in Fig. 7 together with the 
excluded a peak of ^^'^W. 

Despite the obtained limit is far away from the theoretical predictions (f.i., 
Ti/2 ~ 1.3 X 10^° yr [63]), the limit is almost two orders higher than the previous 
one Ti/2 > 1-0 x 10^^ yr derived from the low background measurements with a small 
(4.5 g) ZnW04 crystal scintillator |64j. 

4. Conclusions 

A low background experiment to search for 2/3 processes in ®^Zn, "^^Zn, ^^°W, and ^^^W 
was carried out over more than 19 thousands hours in the underground Gran Sasso 
National Laboratories of the INFN by using radiopure ZnW04 crystal scintillators. The 
total exposure of the experiment is 0.5295 kg x yr. 

New improved half-life limits on double electron capture and electron capture with 
positron emission in ®^Zn have been set in the range: 10^^ yr to 10^^ yr depending on the 
mode. The indication on the {2i>-\-0i>)eP~^ decay of ^^Zn with T1/2 = (1.1±0.9) xlO^^ yr 
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suggested in [58 j is completely disproved by the results of the present experiment. Note 
that to date only four nuclides ^^Ca, ^^Kr, ^^^Sn, and ^^°Te) among 34 candidates to 2e, 
ef3~^, and 2/3"'" processes were studied at similar level of sensitivity in direct experiments. 
However, it is worth noting that the limits are still far from theoretical predictions. 

In addition to ^^Zn decays, in the course of the present experiment two important 
by-products were obtained: (1) the new half-life limits on the 2/3 processes in ^°Zn, 
^^°W, and ^'^^W on the level of 10^^ — 10^^ yr (the 0i'2e capture in ^'^^W is of particular 
interest due to the possibility of resonant process); (2) rare a decay of ^^^W with a half- 
life Ti/2 = 1.31q;5 ^ 10^^ has been observed and new half-life limit on a transition of 
^^^W to the 1/2- 375 keV metastable level of ^^^Hf has been set as T1/2 > 6.7 x 10^° yr. 
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